The Annals of Applied Probability 
2008, Vol. 18, No. 1, 59-99 
DOI: 10.1214/07-AAP444 

© Institute of Mathematical Statistics. 2008 



"oo". VOTER MODELS WITH HETEROZYGOSITY SELECTION 

o . 

O ■ By Anja Sturm 1 and Jan Swart 2 



(N 
c3 



in 

(N 



in 
in 



o 



University of Delaware and UTIA Prague 



This paper studies variations of the usual voter model that favor 
types that are locally less common. Such models are dual to cer- 
tain systems of branching annihilating random walks that are parity 
preserving. For both the voter models and their dual branching an- 
nihilating systems we determine all homogeneous invariant laws, and 
we study convergence to these laws started from other initial laws. 

Oj 

_C 1. Introduction and main results, 

-t— > 

1.1. Voter models with heterozygosity selection. This paper studies vari- 
ations of the usual voter model that favor types that are locally less common. 
These systems can be used to model the distribution of two types of organ- 
isms (two similar species or merely different genetic variants of the same 
*n . species) that occupy overlapping ecological niches, and therefore compete 

with each other for resources. If both types are equally fit, but their ecolog- 
ical niches are not completely identical, then individuals belonging to the 
type that is locally less common have an advantage, since they can use re- 
sources that are not used by most of their neighbors. This effect is called 
negative frequency dependent selection or (positive) heterozygosity selection. 
(Here, following a common practice in population biology, the word het- 
erozygosity refers to the degree of genetic variation within a population as 
a whole, rather than the variation between homologous chromosomes in a 
J> , diploid organism.) 

' Our processes of interest are Markov processes X = (X t )t>o with state 

5_i 1 space {0,l} zd . We denote a typical element of {0, l} zd by x = (x(i)) i6 
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where x{i) € {0, 1} is interpreted as the type of the organism at the site i. 
Borrowing terminology from physics, we sometimes also call x(i) the spin 
at i. We say that a {0, l} zd -valued Markov process X is spin-flip symmetric 
if its dynamics are symmetric under a simultaneous flip of all spins, that is, 
the transition x i— ► x' happens at the same rate as the transition (1 — x) i— > 
(1-x'). 

More specifically, we are interested in the following models. 

Definition 1 {Neutral Neuhauser-Pacala model). The neutral Neuhauser- 
Pacala model is the spin-flip symmetric Markov process in {0,1} Z such 
that if the state of the process is x, then x(i) flips from to 1 with rate 
/i(/o + otfi) (and likewise for jumps from 1 to 0, by spin-flip symmetry), 
where 

(1-1) fr-= jhl E l{.<a=r} <T = 0,1) 

denotes the local frequency of type r in the block of (2R + l) d - 1 sites 
centered around i, not containing i itself, given by 

(1.2) Mi := {j = (ji, . . . ,j d ) £ Z d : < \j k - i k \ < R Vfc}, 

and < a. < 1, d, R > 1 are parameters such that minjii, i?} > 2. 

The neutral Neuhauser-Pacala model is a special case of the model intro- 
duced in [21], when their parameters satisfy A = 1 and qoi = aio =: ot (called 
the "symmetric case" there). For a = 1, this is a usual (range R) voter model, 
while for a < 1, locally rare types have an advantage. Neuhauser and Pacala 
interpret the rate /i(/o + otfi) as follows: at each site i, an organism of 
type dies with rate /o + af\ due to competition with its neighbors, and 
is immediately replaced by an organism of a random type chosen from Mi. 
If a < 1, then the interspecific competition is smaller than the intraspecific 
competition, hence, locally rare types die less frequently. 

Our next model of interest is another nonlinear voter model. 

Definition 2 (Affine voter model). The affine voter model is the spin- 
flip symmetric Markov process in {0, 1} Z such that if the state of the process 
is x, then x(i) flips from to 1 with rate afi + (1 — a)l{j 1> o}, where fo, f\ are 
defined as for the previous model. Here < a < 1 and d,R> 1, min{d, R} > 
2. 

Again, for a = 1, this is a usual (range i?) voter model, while for a < 1, 
locally rare types have an advantage. For a = 0, the affine voter model is 
a threshold voter model [7, 14, 18]. In this and the previous definition, we 
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have excluded the case d = 1 = R since this model has special behavior (see 
Section 2.1). 

The third and final model we will consider is a one-dimensional model. 

Definition 3 (Rebellious voter model). The rebellious voter model is 
the spin-flip symmetric Markov process in {0, 1} Z such that if the state of 
the process is x, then x(i) flips with rate 

a(l{x(i-i)#x(*)} + Mx(i)^x(i+i)}) 

(l.oj 

+ (1 - a)(l{a;(i_2)^x(i-l)} + l {x{i+l)^x(i+2)})- 

For a = 1, this is a nearest-neighbor one-dimensional voter model, but 
for a < 1, locally rare types have an advantage. To see this, note that if 
x(i + 1) ^ x(i + 2) and the spin at i flips, then x(i) was previously of the 
most common type in the set {i,i+ 1, i + 2}. 

In all of the models above, we call a the (interspecific) competition param- 
eter. As usual with voter models, the main interest in these models lies in 
the phase transition between coexistence and noncoexistence. It is believed, 
and has been proved in special cases, that coexistence occurs for sufficiently 
small competition parameters a or in high dimensions, while noncoexistence 
occurs in low dimensions and for large competition parameters. 

We note that Blath, Etheridge and Meredith [1] have studied systems of 
interacting Wright-Fisher diffusions with heterozygosity selection, that is, 
the stepping stone analogue of the voter models discussed in the present 
paper. 

1.2. Duality with parity preserving processes. The neutral Neuhauser- 
Pacala model and the affine and rebellious voter models are cancellative 
spin systems. Here, a Markov process X = (Xt)t>o with state space {0, l} z<i 
is a cancellative spin system if its dynamics are of the following special form. 
For each finite set AcZ d x there is a rate a(A) > 0, such that with this 
rate, the process jumps from the state x to x + Ax mod(2), where 

(1.4) Ax{i):= x{j) mod(2) (ieZ d ). 

r-(i,j)eA 

(In [13], a somewhat more general class of cancellative spin systems is con- 
sidered, where also spontaneous flips are allowed.) 

The neutral Neuhauser-Pacala model can be cast into the general form 
of a cancellative spin system, with rates a(A) given by 

ct 

(1-5) [_ a 

a({i} x {j,k}) = (j,k£Afi, j ^ k) 
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and a{A) = in all other cases. Likewise, the affine voter model can be 
formulated as a cancellative spin system, with rates a{A) given by 

a{{i} x {i,j}) =a\M l \- 1 (je^Vi), 

(1.6) 

a({i} xA) = (l-a)2-l Af 'l +1 (A C M U {i}, |A| even), 

while in case of the rebellious voter model, the rates are given by 

. , a({i} x {i-l,i}) = a({i} x {i,i + l}) = a, . . 

{ '> a({i}x{i-2,i-l}) = a({i}x{i + l,i + 2}) = l-a, lt ' 

and a{A) = in all other cases. 

Under a suitable summability assumption on the rates, for each cancella- 
tive spin system X there exists a unique cancellative spin system Y such 
that X and Y are dual to each other in the sense that 

(1.8) mX t Y \ is odd] = F[\X Y t \ is odd] (t > 0), 

whenever X and Y are independent (with arbitrary initial laws), and either 
|Xo| or |Yo| is finite (see [13]). Here, for any x,y € {0, 1} Z , we write \x\ := 
Y^i x (i) an d (xy)(i) := x(i)y(i). If X is defined by rates ax {A), then Y is de- 
fined by rates ay {A) given by ay {A) = ax(A T ), where A T := : (i, j) € 
A}. We note that since the functions y i— > li\ xy \ i s odd} with |x| < oo are dis- 
tribution determining, (1.8) determines the transition laws of Y uniquely. 

For example, if X is the rebellious voter model, then the dynamics of Y 
have the following description. We interpret the sites i for which Yt(i) = 1 
as being occupied by a particle at time t. Then each particle jumps with 
rate a one step to the left, and with the same rate to the right, with the 
rule that if the site the particle lands on is already occupied, the two parti- 
cles annihilate. Moreover, with rate 1 — a, each particle gives birth to two 
new particles located on its nearest and next-nearest site to the left, and 
with the same rate on the right, again annihilating with any particles that 
may already be present on these sites. We call Y the asymmetric double 
branching annihilation random walk (ADBARW). Also in case of the neu- 
tral Neuhauser-Pacala and affine voter models, one may check that the dual 
model is a system of branching annihilation random walks, with branching 
rate proportional to 1 — a. Since the number of particles that are born is 
always even, these systems preserve parity. 

More generally, if A is a cancellative spin system and Y is its dual, let 
us say that Y is parity preserving if the process started in any finite initial 
state satisfies (— l)' y '' = (—1)'^°' a.s. It is not hard to see that the follow- 
ing statements are equivalent: (i) X is spin-flip symmetric, (ii) Y is parity 
preserving, (iii) ax (A) = unless \{j : (i,j) £ A}\ is even for all i G 

We next define the concepts of coexistence, persistence, survival and stabil- 
ity. Below, and 1 denote the configurations in {0, l} zd which are identically 
zero and one, respectively. 
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Definition 4 (Coexistence). We say that a probability law \i on {0, 1} 
is coexisting if /i({0, 1}) = 0. We say that a spin-flip symmetric cancellative 
spin system X exhibits coexistence if there exists a coexisting invariant law 
for X. 

Definition 5 (Persistence). We say that a probability law \i on {0, l} zd 
is nonzero if /x({0}) = 0. We say that a parity preserving cancellative spin 
system Y exhibits persistence if there exists a nonzero invariant law for Y. 

Below, X x and Y v denote the processes X and Y started in Xq = x and 
Yq = y, respectively. 

Definition 6 (Survival). We say that a spin-flip symmetric cancellative 
spin system X survives if 

(1.9) P[Xf ^ V< > 0] > for some \x\ < oo. 

We say that a parity preserving cancellative spin system Y survives if 

(1.10) P[Y t y ^ Vt > 0] > for some |y| < oo, \y\ even. 

Note that because of parity preservation, the left-hand side of (1.10) is 
always one if \y\ is odd. In view of this, it will usually be clear what we 
mean when we say that a process survives. Our definitions are not entirely 
unambiguous, however, since it is possible for a cancellative spin system to 
be both spin-flip symmetric and parity preserving. When there is danger of 
confusion, we will say that the even process survives if we mean survival in 
the sense of (1.10). 

For a parity preserving cancellative spin system Y started in an odd initial 
state, we define 

Y t (i):=Y t (l(t)+i) (ien d ) 

(1-11) 

where l(t) := inf{z € Z d : Y t (i) = 1}. 

Here, the infimum is defined componentwise. Note that {Yt)t>Q is a Markov 
process with state space {y € {0, l} Nd : \y\ is finite and odd, inf {i : y(i) = 1} = 
0}. We call Y the process Y viewed from its lower left corner. 

Definition 7 (Stability). We say that a parity preserving cancellative 
spin system Y is stable if the state with one particle at the origin is positively 
recurrent for the Markov process Y . 
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Using duality, one can prove that each cancellative spin system has a 
special invariant law, which is the limit law of the process started in product 
law with intensity 1/2 [13]. Thus, 



where an d za/ are invariant laws for the processes X and Y, respec- 

1/2 

tively. Because of certain analogies with additive spin systems, we call 

1/2 

and Uy the odd upper invariant laws of X and Y , respectively. 
The next lemma is a simple consequence of duality. 

Lemma 1 (Invariant laws and survival). Let X be a spin- flip symmetric 
cancellative spin system and let Y be its dual parity preserving cancellative 
spin system. Then: 

(a) The following statements are equivalent: (i) X exhibits coexistence, 

1/2 

(ii) is not concentrated on {0,1}, (hi) Y survives. 

1/2 

(b) The following statements are equivalent: (i) Y persists, (ii) Vy is 
not concentrated on 0, (hi) X survives. 

Apart from the relations in Lemma 1, for our models of interest, one 
readily conjectures a number of other relations between the concepts we 
have just introduced. While being supported by numerical simulations, these 
conjectures appear to be hard to prove in general. Thus, we formulate as 
open questions: 

Ql. Is coexistence of X equivalent to survival of XI 

Q2. Does stability of Y imply extinction of Y? 

Q3. Does coexistence for a imply coexistence for all a' < a? 

The rebellious voter model has a special property, explained in Section 2.1, 
that allows us to answer Question Ql positively for this model. 

Lemma 2 (Survival and coexistence). The rebellious voter model sur- 
vives if and only if it exhibits coexistence. 

1.3. Results. We say that a probability law fi on {0, l} zd is homogeneous 
if fj, is translation invariant. The odd upper invariant laws v^ 2 and Vy 2 are 
examples of homogeneous invariant laws, and so are the delta-measures do 
and (in case of X) Si- We will see that, under weak additional assumptions, 
there are no others. 




(1.12) 
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For additive spin systems, duality may be employed to show that, under 
certain conditions, the upper invariant law is the only nonzero homogeneous 
invariant law, and the long-time limit law for any system started from a 
nonzero homogeneous initial law; see [15] and [17], Theorem III. 5. 18. With 
certain complications, these techniques can be adapted to cancellative spin 
systems and their odd upper invariant laws. This idea has been success- 
fully applied in [4] to certain annihilating branching processes. As our next 
theorem demonstrates, it can be made to work for our models as well. 

Theorem 3 (Homogeneous invariant laws). Let X be either the neutral 
Neuhauser-Pacala model, the affine voter model or the rebellious voter model 
and let Y be its dual. Then: 

1/2 

(a) If a < 1 and Y survives, then u-J is the unique homogeneous coex- 
isting invariant law of X. If, moreover, a > and Y is not stable, then the 
process X started in any homogeneous coexisting initial law satisfies 



(b) If a < 1 and X survives, then u Y is the unique homogeneous nonzero 
invariant law ofY. If, moreover, a > and d>2, then the process Y started 
in any homogeneous nonzero initial law satisfies 



Theorem 3 does not tell us anything about the values of a for which 
coexistence occurs. For the affine voter model, coexistence for a = has 
been proved in [7, 18]. It seems likely that, using comparison with oriented 
percolation, this result can be extended to small positive a. Comparison 
with oriented percolation was used in [21], where coexistence for the neutral 
Neuhauser-Pacala model for a sufficiently close to zero is proved, and in [9] 
where (among other things) coexistence is proved for the neutral Neuhauser- 
Pacala model in dimensions d > 3 for a sufficiently close to one. 

Again using comparison with oriented percolation, we will prove that, for 
small a, the rebellious voter model coexists. In fact, we prove considerably 
more. 

Theorem 4 (Complete convergence). Let X be the rebellious voter model 
and let Y be its dual. Then there exists an a' > 0, such that for all a 6 



(1.13) 




(1.14) 




t— »oo 



[0,a'): 



(a) The process X exhibits coexistence and survival. The process Y ex- 
hibits persistence, survival and is not stable. 
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(b) The process X started in an arbitrary initial law satisfies 
(1.15) F[X t G •] Mo + pxffj + (1 - po - pi)4 /2 , 

r — >oo — — 

where p q := P[X t = q for some t>0] (q = 0, 1). 

In analogy with similar terminology for the contact process, we call the re- 
sult in part (b) complete convergence. Complete convergence for the thresh- 
old voter model (i.e., the affine voter model with a = 0) was proved in [14]. 

2. Methods and discussion. 

2.1. Interfaces. For one-dimensional models, there is, apart from duality, 
a useful additional tool available. If X is a spin-flip symmetric cancellative 
spin system on Z, then setting 

(2.1) Y t {i) := l {Xt(i , A ( i+1)} (t > 0, i G Z) 

defines a Markov process y = (Yt)t>o in {0, 1} Z that we call the interface 
model associated with X. Under a suitable summability assumption on the 
rates, Y is a parity preserving cancellative spin system. 

By looking at interfaces, we can explain our interest in the rebellious voter 
model. Moreover, we can explain why we have excluded the case d = 1 = R 
from Definitions 1 and 2. Consider the following spin-flip symmetric models. 



Definition 8 {One- dimensional models). For any < a < 1, the dis- 
agreement voter model is the cancellative spin system on Z with rates a(A) 
given by 

(2 2 \ °( W * {» - 1» *}) = a ( W x {h i + 1}) = a, , _ 

^ ' ' a({i}x{i — l,i + l}) = l — a, ^ 

and the swapping voter model is given by the rates 

(2 31 «({«} x - 1,^}) = a({i} x + 1}) =a, 

^ ' ' a({i,i + 1} x + 1}) = 1 - a, ^ '' 

In each case, it is understood a(A) = for all A other than those mentioned. 

If in the definition of the neutral Neuhauser-Pacala model with competi- 
tion parameter onp one would set d = 1 = R, then up to a trivial redefinition 
of the speed of time, one would obtain a disagreement voter model with pa- 
rameter a = 2«np/(1 + oinp)- Likewise, setting d = 1 = R in the definition 
of the affine voter model with competition parameter oav yields, up to a 
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change of the speed of time, a disagreement voter model with parameter 
a = 1/(2 - am)- 

To explain the special properties of the disagreement and rebellious voter 
models, we look at the way these models are related to other models through 
duality and interface relations. These relations are summarized in Figure 1. 
Recall that the dual of the rebellious voter model is the ADBARW. The 
dual of the disagreement voter model has been called the double branching 
annihilating random walk (DBARW) in [24]. The swapping voter model has 
a mixture of voter model and exclusion process dynamics. We call its dual 
the swapping annihilating random walk (SARW). 

The swapping annihilating random walk (SARW) has the special prop- 
erty that the number of particles cannot increase. As a result, the behavior 
of the disagreement and swapping voter models and their duals is largely 
known. For each a > 0, the disagreement and swapping voter models exhibit 
extinction and noncoexistence. The DBARW and the SARW get extinct and 
are not persistent. These facts follow from [24], Theorem 8 and [21], Theo- 
rem 2(b). The disagreement voter model with a = 1/2 is a threshold voter 
model and has been studied earlier in [7]. It is trivial that the SARW is 
stable. We will prove stability for the DBARW elsewhere [23]. 

The rebellious voter model has the special property that its interface 
model and dual coincide, which is one of our main reasons for introducing 
it. This property is very helpful in the proof of Theorem 4. Moreover, it 
allows us to prove Lemma 2: 

Proof of Lemma 2. It is easy to see that X survives if and only if 
its (even) interface model Y survives. Since Y is also the dual of X, by 
Lemma 1, X exhibits coexistence if and only if Y survives. □ 

For the rebellious voter model, numerical simulations suggest the exis- 
tence of a critical value a c ~ 0.5 such that, for a < a c , one has survival, 
coexistence and instability of the interface model, while, for a > a c , one 
has extinction, noncoexistence and stability of the interface model (see Fig- 
ure 2). We conjecture qualitatively similar behavior (but with a different 
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Fig. 1. Relations between models. 
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Fig. 2. Equilibrium density of the ADBARW as a function of a. The data were obtained 
starting with one particle on an interval of 700 sites with periodic boundary conditions, 
gradually lowering a from, one to zero during a time interval of length 300, 000. 

critical point) for the one-dimensional neutral Neuhauser-Pacala and affine 
voter models. Note that in this respect, these models are quite different from 
the disagreement voter model, which gets extinct for all a > 0. 

2.2. Homogeneous laws. In this section we discuss the methods used to 
prove Theorem 3. Using duality, it is not hard to see (see Section 3.1) that 
the odd upper invariant law of X, defined in (1.12), is uniquely determined 
by 

(2.4) P[|^ /2 y| isodd] = iP[r/#0Vs>0] (|l/|<oo). 

In order to prove the convergence in (1.13), it therefore suffices to show that 

(2.5) lim F[\X t y\ is odd] = ±P[Y/^0 Vs > 0] (\y\<oo). 

t — ►oo 

Using duality (1.8), we can rewrite the left-hand side as 

(2.6) lim n\Xt Y t \ | is odd], 

where to > is fixed. Conditioning on the event of survival, we need to show 
that 

(2.7) lim nW t y _ t0 | is odd 1 17 ^ Va > 0] = \. 

It turns out that we can show this, provided that we can show that 

(2.8) lim P[0< \Y?\ < N] = (N>1), 

t— >oo 

that is, Y exhibits a form of extinction versus unbounded growth. To prove 
this, we need to assume that Y is not stable. The proof uses induction on 
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N and is quite subtle. Originally, we only knew how to prove convergence 
in Cesaro mean, until we saw the paper of Handjani [14] where usual con- 
vergence as t — > oo is proved in the context of threshold voter models. 

The proof of the convergence in (1.14) follows similar lines. This time, 
instead of (2.7), we need to prove 

(2.9) lim n\X?_ t0 Y t0 \ is odd | X* + Vs > 0] = ±. 

At first, it might seem that this is true provided that, in analogy with (2.8), 
one has that P[0 < \Xf\ < N] -> for all N > 1. Indeed, the reference [22] 
contains a claim of this sort, but as we explain in Section 3.6 below, this 
is not correct. Indeed, in order for the probability on the left-hand side of 

(2.9) to be close to ^, we need many events that could affect the parity of 
|^f-t ^ol- This means that there must be many sites i,j,k, close to each 
other, such that Yq(z) = 1, while Xf_ tQ (j) ^ Xf_ tQ {k). Indeed, it suffices to 
verify that, conditional on survival, the quantity 

(2.10) |VAf | := eZ d xZ d :\i-j\= l,Xf(i) + Xf(j)}\ 

tends to infinity as t — > oo. In dimensions d>2, we can verify this, but in 
dimension d = 1, we run into the difficulty that it is hard to rule out the 
scenario that, at certain large times, Xf consists of just one large interval 
of ones, in which case |VAf | = 2. (Indeed, to prove a result in d = 1, one 
would probably need to assume that the interface model associated with X 
is not stable. We have not carried out this approach.) 

2.3. Comparison with oriented percolation. Theorem 4 is similar to the 
main result of [14]. In that paper a key technical tool is comparison of 
threshold voter models with threshold contact processes, which in [18] were 
shown to survive. Using complete convergence for these threshold contact 
processes, one can then prove the analogue statement for the threshold con- 
tact processes. 

Our approach will be similar, except that we will use comparison with 
oriented percolation. Also, instead of proving a comparison result for the 
rebellious voter model X , we prefer to work with the ADBARW Y, which 
is both the dual and the interface model associated with X . We will prove 
that if a is small, then for each p < 1, the process Y viewed on suitable 
length and time scales dominates an oriented percolation with parameter p. 
This kind of rescaling argument was first used in Bramson and Durrett [3] . 
In [4, 11], it was shown how the technique may be amended to cover also 
interacting particle systems that lack monotonicity. Unfortunately, these 
references are somewhat imprecise when it comes to showing m-dependence 
for the objects one compares with; this is done more carefully in [12]. We 
will use a somewhat different argument to ensure m-dependence than the 
one used in that reference. 
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Since the comparison result is of some interest on its own, we formulate 
it here. We first introduce oriented (site) percolation with percolation pa- 
rameter p. Let Zgygu := {(x, n) E Z 2 : x + n is even}. Let {u> z : z E Z 2 vcn } be 
i.i.d. Bernoulli random variables with F[lu z = 1] = p. For z,z' E Z 2 vcn , we 
say that there is an open path from z to z' , denoted as z — ► z', if there 
exist (x n ,n), . . . , (x m ,m) E Z 2 vcn with \xk — Xk-i \ = 1 and U( Xk ,k) = 1 f° r all 
n < k < m, such that (x n ,n) = z and (x m ,m) = z' . By definition, z — > z for 
all z E ^even- For given ^4 C Z even := {2n : n E Z}, we put for n > 

(2.11) W n := {x E Z : (x, n) E Z 2 vcn , 3x' E A s.t. (a/, 0) -> (x, n)}. 

Then = (W n ) n >o is a Markov chain, taking values, in turn, in the subsets 
of Z even and Z 0< jd, started in Wq = A. We call W the oriented percolation 
process. 

The comparison entails defining certain "good" events concerning the 
behavior of the ADBARW in large space-time boxes. Let L > 1 and T > 0. 
For any x E Z, put 

J a . := {2Lx — L, . . . , 2Lx + L} and 

(2.12) 

4 := {2Lx - 4L, . . . , 2Lx + 4L}. 

Let Y be an ADBARW started in an arbitrary initial state Yq. We define a 
set of "good" points for n > by 

Xn ■= {x E Z : (x, n) E Z^ ven , 3i E 4 s.t. y„T(«) = 1 and, 

(2.13) 

in case n > 1, V(n -1)T <t< nT 3i E 4 s.t. F t (i) = 1}. 
With these definitions, our result reads: 

Theorem 5 (Comparison with oriented percolation). For each p < 1, 
there exists an a' > suc/i £/ta£ /or all a E [0, a') there exist L > 1 and 
T > 0, stic/i i/iai «/ Y is an ADBARW with parameter a, started in an 
arbitrary initial state Yq, then the process (xn)n>o defined in (2.13) can be 
coupled to an oriented percolation process (W n ) n >o with parameter p and 
initial state Wq = xo, in such a way that W n C Xn for all n > 0. 

2.4. Discussion and open problems. Three open problems have already 
been formulated in Section 1.2. As a fourth problem, we mention the follow- 
ing: 

Q4. For which values of the parameters do voter models with heterozygosity 
selection exhibit noncoexistence? 
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As mentioned before, there are several results proving coexistence for the 
sort of models we are considering. Very little is known about noncoexistence, 
except for the pure voter model in dimensions d = 1, 2 and the disagreement 
voter model, which is somehow special. There is some hope that the methods 
used in [9] (see also [8]), who prove coexistence for the neutral Neuhauser- 
Pacala model with a close to one in dimensions d > 3, can be extended to 
cover dimension 2 as well. Therefore, noncoexistence can be expected in di- 
mension one only. In fact, it seems that physicists believe that for these and 
similar models, there is a critical dimension d c ~ 4/3 such that only below 
d c there is a nontrivial phase transition between coexistence and noncoexis- 
tence [5, 6, 25]. Here, the fractional dimension probably refers to self-similar 
lattices. 

Note that, by duality, proving noncoexistence for large interspecific com- 
petition boils down to proving extincton for a parity preserving branching 
process with a small branching rate. Usually, for interacting particle systems, 
it is easier to find sufficient conditions for extinction than for survival; this is 
true for the contact process, and also for the annihilating branching process 
studied in [4]. The difficulties in our case come from parity preservation, 
which makes extinction difficult, slow (slower than exponential) and "nonlo- 
cal," since it may require particles to come from far away to annihilate each 
other. 

2.5. Outline. The rest of the paper is devoted to proofs. Lemma 2 has 
already been proved in Section 2.1 above. In Section 3 below, we prove 
Lemma 1 and Theorem 3. Theorems 4 and 5 are proved in Section 4. 

3. Homogeneous invariant laws. 

3.1. Generalities. In order to prepare for the proof of Theorem 3, we 
will derive some results for a general class of cancellative spin systems on 
lj d . Throughout this section, X is a cancellative spin system on Z rf defined 
by rates a(A) =ax(A) (see Section 1.2), and Y is its dual cancellative spin 
system defined by rates ay (A) = ax(A T ). In order for X and Y to be well- 
defined, we make the summability assumptions: 



We also assume that our rates are (spatially) homogeneous, in the sense that 





sup$^a(A)|{i:(i,j) € A}\ < oo, 

1 A3i 

sup^a(A)|{j:(j,i) G A}\ < oo. 



1 ABi 
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We will sometimes need the graphical representation of X. Independently 
for each finite A 6 Z rf x let tta be a random, locally finite subset of M, 
generated by a Poisson processes with intensities a{A). We visualize this by 
drawing an arrow from (i,t) to (j,t) for each t 6 71^4 and 6 ^4 (note the 
order). By definition, a pat/j from a subset C C Z rf x R to another such subset 
D is a sequence of points iq, . . . , i n S Z rf and times to < *i < • • ■ < t n < t n +i 
(n > 0) with 

(3.3) VI < m < n 3A s.t. t m 6 7r^ and (i m ,i m _i) € .A, 

such that (iojio) € C and (i n ,t n+ i) G -D. Thus, a path must walk upward 
in time and may jump from one site to another along arrows. With these 
conventions, for any subset U C Z rf , setting for t > 0, i € Z d , 

(3.4) -^t(*) := l{thc number of paths from Ux{0} to (i,t) is odd} 

defines a version of the cancellative spin system X defined by the rates a(X), 
with initial state Xq = Ijj. 

Proof of Lemma 1. We only prove part (a), the proof of part (b) 
being similar. By duality (1.8), for each y £ {0, l} zd with \y\ < 00, 

(3.5) n\xl /2 y\ is odd] = P[|X 1/2 y/| is odd] = ±P[lf + 0], 
hence, 

(3.6) mx\ l2 y\ is odd] — ► ±P[Y/ + Vt > 01. 

Recall that the functions x 1— > l/u^i i s odd} with |y| < 00 are distribution 

determining. Therefore, the odd upper invariant law of X is uniquely deter- 

1/2 

mined by (2.4). If Y survives, then (2.4) shows that W[\XoL y\ is odd] > 

for some \y\ even, hence, za!/ 2 is not concentrated on {0,1}. This shows 

1/2 

that (iii)=>(ii). To show that (ii)=>(i), it suffices to note that P[Xoo £ • | 
1/2 

^00 7^0>i] is a coexisting invariant law. To see that (i)=^(iii), assume that 
P[^oo € •] is a coexisting invariant law for X. Let Si 6 {0, l}^ d be defined by 
5i(j) = 1 if i = j and otherwise. Then by the duality (1.8) applied to the 
process X started in Xoo, we have, for i ^ j,t > 0, 

(3.7) F[Y* i+5j + 0] > FUX^Yf^ I is odd] = P[X oc (i) + X^j)} > 0, 
and therefore, 

(3.8) P[y/ l+5j ± Vt > 0] = lim P[y/ l+<5j ^ 0] > P[X oc (i) + X^j)} > 0, 
which shows that Y survives. □ 
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For later use, we introduce some more notation. The set {0, 1}, equipped 
with multiplication and addition modulo 2, is a finite field. We can view 
{0, 1} Z as a linear space over this field. In this point of view, if we identify 
a finite set A C l> d x Z d with the matrix A such that A(i,j) = 1 if G A 
and A(i,j) = otherwise, then Ax as defined in (1.4) is the usual action of 
a matrix on a vector: 

(3.9) Ax(i)=^2A(iJ)xU) mod(2). 

3 

In line with these observations, for x G {0, l} zd , we let x T denote the [mod(2)] 
linear form on {0, l} z<i given by 

(3.10) x T y :=^2x(i)y(i) mod(2) = l^ is odd} , 

i 

which is well defined whenever \xy\ < oo. 

3.2. Uniqueness and convergence. In this section we continue to work 
in the general set-up introduced above. We show how a sort of "extinction 
versus unbounded growth" for the dual process Y can be used to prove 
convergence to equilibrium for X , started in a homogeneous initial law. Our 
result is similar in spirit to claims by Simonelli [22] . Unfortunately, as already 
mentioned in Section 2.2, that reference contains an error, which we point 
out in Section 3.6 below. 

To formulate our result, we need to identify "good" configurations where 
parity can change. To this aim, we select a finite set B whose elements are 
finite subsets B of Z rf x Z d such that a(B) > 0, and we define, for y G {0, 1} Z , 

(3.11) \\y\\ B :=\{ieZ d :3xe {0, 1} Z< \ B G B s.t. y T {TiB)x = 1}\. 

Note that y T (TiB)x = 1 is equivalent to \{(j,k) G TiB:y(j) = 1 = x(fc)}| 
being odd, or, equivalently, to \x'y \ having a different parity from \xy\, where 
x' = x + (TiB)x mod(2). 

For our models of interest, we may choose || • \\q as follows. If X is a neutral 
Neuhauser-Pacala model affine voter model or rebellious voter model (with 
arbitrary a), then we can find i,j such that 

(3.12) a({0}x{i,j})>0. 

Taking for B, the one-point set {{0} x {i,j}} now leads to ||y||s = \y\- If, 
on the other hand, X is the dual of any of these models, then for each 
k = 1, . . . , d, we can find such that 



(3.13) 



a({0,e fc } x {i k })>0, 
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where e\ := (1, 0, . . . , 0), . . . , := (0, . . . , 0, 1) are unit vectors in each of the 
d dimensions of 7L d . Taking for £>, the set {{0, e\} x {i\}, ■ ■ ■ , {0, e^} x {id}} 
now yields = \Vy\ [recall (2.10)]. 

We say that a probability law fi on {0, 1} Z is locally nonsingular if, for 
every finite AcZ rf and every j/G{0,l} A , 

(3.14) n({x:x{i)=y{i) Vi G A}) >0, 

that is, every finite configuration has positive probability. We say that ^ is 
X-nontrivial if \x is concentrated on states x such that PLYf G •] is locally 
nonsingular for each t > 0. 

Theorem 6 (Homogeneous invariant laws). Assume that X has a ho- 
mogeneous X-nontrivial invariant law. If the dual process Y v started in any 
finite initial state YJf = y satisfies 

(3.15) F[3t >0 s.t. \\Y t y \\ B i {l,...,iV}] = 1 (N>1, \y\ <oo), 

then z^ x is X-nontrivial, and the unique homogeneous X-nontrivial invari- 
ant law of X . If, moreover, one has 

(3.16) limP[||y i y || B ^{l,...,iV}] = l (JV>l,M<oo), 

t— >oo 

then the process X started in any homogeneous X-nontrivial initial law sat- 
isfies 

(3.17) F[X t G •] => 

t — ^oo 

The proof depends on two lemmas. 

Lemma 7 (Parity uncertainty). For each e > and t > 0, there exists 
an N > 1 such that if X and Y are started in deterministic initial states Xq 
and Yq satisfying 

(3.18) \{i £Z d :3B £B s.t. Y^(TiB)X = 1}\ > N, 
then 

(3.19) |P[|-Xt*ib| ^ odd] - \\ <e. 

Proof. Let A denote the set of all finite subsets of Z d x Z d . For any 
A G A, let us put 

R_(A) := {j : A(i,j) = 1 for some i G 

(3.20) 

R+(A) := {i : A(i, j) = 1 for some j G Z d }, 
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and let us call R{A) := R_{A) U R+(A) the range of A. Let us say that 
A, A G A are disjoint if i?(^4) n R(A) = 0. If X ,Y satisfy (3.18), then 
it is not hard to see that we can successively choose i±,...,i n G 7L d and 
Bi,...,B n &B such that T^Bi, ...,T in B n are disjoint, Y^(T im B m )X = 1 
for all 1 < m < n, and 

(3.21) n>N/(\B\K 2 ), 

where K := m<%x{\R(B)\ : B G B}. Indeed, for each B,B' eB and i G Z d , 
there are at most K 2 points j G Z d such that R(TjB') n R&B) ^ 0, so 
once we pick a point from the set in (3.18), there are at most \B\K 2 points 
we cannot pick anymore. We now use the graphical representation (see Sec- 
tion 3.1). Let 

M := {m : 1 < m < n, W im B m D (0, t)\ G {0, 1} and 

(3.22) \ir A n (0,t)\ =0VAeA,A^ T im B m 

s.t. R(A)DR(T im B m )^0}. 

Thus, for each m G M, the Poisson process associated with Ti m B m becomes 
active zero or one time during the time interval (0, t), and no other Poisson 
process creates arrows in the range R(Ti m B m ) during this time interval. 
We claim that if N is sufficiently large, then the set M is large with high 
probability. 

Indeed, since a(Ti m B m ) = a(B m ) is bounded from above, the probability 
that \iTT im B m H (0, <)| G {0, 1} is uniformly bounded from below. By summa- 
bility [condition (3.1)], the probability that no other Poisson process creates 
arrows in the range R(Ti m B m ) during (0, t) is also uniformly bounded from 
below. These events are not independent for different m, but they are pos- 
itively correlated, so we get a lower bound assuming independence, which 
proves our claim. 

Let T be the cr-field generated by the random set M and by all Poisson 
processes tta H (0, t) with A G A\{Ti m B m :m G M}. Thus, T corresponds 
to knowing the random set M and all Poisson processes on (0,i), except 
those associated with the Ti m B m with m G M. Set m := a(B m )t. Note that 
the 9 m are uniformly bounded from above and below by the fact that B is 
finite. We claim that if we condition on then under the conditioned law, 
the random variables \iTTi B m H (0,t)| are independent {0, l}-valued random 
variables with 



(3.23) P[|vr TimBm n(0,t)| = l|^] = 




(m G M). 



Indeed, first condition on all tta H (0, t) with A ^ T^Bi, . . . ,Ti n B n . Under 
this conditional law, the ttt^Bi H (0, t), . . . ,TTT in B„ H (0,t) are independent 
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Poisson processes with intensities a{B\), . . . , a(B n ). Let 

M' :={m:l<m<n, \tt a n (0,t)| = G ^{T^Bi, . . . ,T in B n } 

(3.24) 

s.t. #(A)ni2(T im S m )/0}. 

Under the conditional law we are considering, M' is a deterministic set, and 
the B m H (0, t) with m 6 M' are independent Poisson processes. Hence, 
since 

(3.25) M = {m G M' : \ir Tim B m D (0, i)| G {0, 1}}, 

if we condition also on the wt- B m H (0,t) with m G M'\M, then under 
this new conditional law, the tt^ B m H (0, £) with m G M are independent 
Poisson processes conditioned to produce zero or one point. This explains 
(3.23). 

Using the graphical representation, we now write Xq = 1u, Yq = ly, and 

(3.26) P[\X t Y \ is even] -F[\X t Y \ is odd] = E[(-l) p+p '], 

where P and P' are the number of paths from U x {0} to V x {t} that 
do and do not use, respectively, arrows created by the Poisson processes 
TtTi B m CI (0, t) with m G M. Note that due to the definition of M the paths 
counted by P use exactly one arrow on (0, t) created by a Poisson processes 
7r T lm B m H (0, i) with m G M. Since Yq - (Ti m B m )Xo = 1 for m G M, this implies 
that P = J2meM Kli B m H (0, t)| mod(2). It therefore follows that 

E[(-l) p+p ' | JP] = (-1) P 'E[(-1) P I F] 

(3.27) =("1) P ' II E[(-l) kT ^ Smn{0,i)l |^] 

= ("1) P ' II ((l-An)-0m), 
mSM 

where the (p m are defined in (3.23). Integrating over the cr-field J-, it follows 
that, under the unconditional law, 



(3.28) |P[|X t y | is even] - F[\X t Y \ is odd]| < E 



n I 1 "2^ 

,m£M 



Since the ^ m are bounded away from zero and one and since \M\ is with 
high probability large if N is large, it follows that for each e > we can 
choose N large enough such that (3.19) holds. □ 

Lemma 8 (Many good configurations). Let X be started in a homo- 
geneous X -nontrivial initial law and let t > 0. Assume that y n G {0, l} z<i 
satisfy 

(3.29) lim \\yn\\B = oo. 
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Then 
(3.30) 



\{ieZ d :3BeB s.t yl{T l B)X t = l}\ 



oo, 



where — > denotes convergence in probability. 



Proof. This proof very closely follows ideas from [15], Theorem (9.2). 



Set 



(3.31) 



C n := {i£Z d :3x£ {0, 1} Z , B £ B s.t. yl{TiB)x = 1}, 
C' n :={ieZ d :3BeB s.t. y^(TiB)X t = 1}. 



By (3.29), \C n 
satisfy 

(3.32) 



oo. We need to show that the random subsets CL C C„ 



s [|c;i>iv]^ 



i 



for all N > 1. By dividing C n into N disjoint sets, each with size tending 
to infinity, we can reduce this to showing (3.32) for N = 1. For each i G C n , 
choose Bi tn G B and Xi jTl G {0, l} zd such that 

(3.33) y£(T i B itn )x i)n = l. 

Since B is finite, by going to a subsequence if necessary, we can assume 
without loss of generality that for some £ G £ and z G {0, 1} R +( B \ where 
R + (B) is as in (3.20), \C n \ — » oo, where 

(3.34) C n := {i G C n :B iin = B, (x i>n (i + j))jeii + (B) = z }- 
It now suffices to prove that P[|C^| > 1] — > 1, where 

(3.35) (5; := {i G C n : (X t (i + j)) jeR+(B) = z}. 
Equivalently, we need to show that 



(3.36) 



E 



IP 



L i&C n 



0. 



Fix e > and k > 1. For each L > 1, we can find, for n sufficiently large, 
subsets C^ k cC n such that > L for all i,j G C^' fe , « / j and | C^' fc | = . 
We claim that there exists an L > 1 such that the process X x started in any 
deterministic initial state X§ = i€ {0, 1} Z satisfies 



E 



{(^(i+j)) jeR+( B)^} 



(3.37) 



n 1 



< II Whmw 



0) 



j6fi + (S) 



7^ 



}] +eA;. 
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One way to see this is to check that the conditions of [17], Theorem 1.4.6 
are fulfilled. Alternatively, one can use the graphical representation. It is not 
hard to see that, for L sufficiently large, the probability that there exist two 
paths between time zero and time t, one ending at R+(TiB) and the other 
at R+(TjB) for some i,j € C~^' k , i ^ j, and both starting at the same site, 
is bounded by ek. This implies (3.37). 

By (3.37), Holder's inequality and the fact that P[Xo £ •] is homogeneous, 
it follows that 



lim supE 

n^oo 

< lim sup f F[X £dx]E 

n— »oo J 



n Hixf 



(i+j))je« + (B)7^} 



(3.38) <limsup/p[X £cte] J] ^{(Xf (B) **}] + ^ 



< 



hmsup [] ( fnX Edx]F[(Xf(i + j)) R 



(B) 



l/k 



+ ek 



= J F[X € dx]F[(X?(j)) jeR+(B) ^ z] k + ek. 

Letting first e — > and then k — > oo, using nontriviality, we arrive at (3.36). 
□ 



Lemmas 7 and 8 combine to give the following corollary. 



Corollary 9 (Parity indeterminacy). Let X be started in a homoge- 
neous X -nontrivial initial law and let t > 0. Assume that y n £ {0, 1} Z satisfy 



(3.39) 
Then 
(3.40) 



lim ||?/n||,6 = °o. 

n— >oo 



lim P[|Xty n | is oc 



Proof of Theorem 6. Imagine that v is a homogeneous X-nontrivial 
invariant law of X. Put 



(3.41) f(y) := / v{dx)l { \ xy \ is odd} (y € {0, 1}* , \y\ < oo). 
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We claim that (f(Y t y )) t >o is a martingale. Indeed, if X is a stationary process 
with law F[X t G •] =v (t GR), independent of Y y , then, by duality (1.8) for 
y G {0, l} zd with |y| < oo, 

(3.42) E[f(Y y )} = F[\X Y y \ is odd] = F[\X t Y y \ is odd] = f(y). 

Using, moreover, the Markov property of Y, this shows that (f (Y y ))t>o is 
a bounded martingale. Set 

(3.43) r N := wf{t > : ||F/|| B £ {1, • . • , iV}}, 

which is a.s. finite for all N > 1 by our assumption (3.15). Hence, by optional 
stopping, 

(3.44) f(y) = E[f(Y y N )] = F[\X t Y y N | is odd| \\Y y N \\ B > iV]P[||Y* || B > N]. 
Letting N — > oo, using Corollary 9, we find that for y G {0, l} zd with |y| < oo, 

(3.45) | is odd} = f(y) = |PKVQ v* > o]. 

1/2 

By (2.4), this implies that = z/ x ' . 

If, moreover, (3.16) holds and X is started in any homogeneous X-non- 
trivial initial law, then, by duality and Corollary 9, 

(3.46) ¥[\X t y\ is odd] = V[\X a Y?_.\ is odd] — ► ^FlY? ^ Vt > 0]. 

1/2 

Since this holds for any finite y, it follows that PpQ G •] => ^ x as t — > oo. 
□ 

3.3. Extinction versus unbounded growth. In order for Theorem 6 to be 
applicable, we need to check that the dual Y of a cancellative spin system X 
satisfies a version of "extinction versus unbounded growth," that is, we must 
show that ||e is either zero or large at random times t or at large fixed 
t. Therefore, in this section, we derive sufficient conditions for a cancellative 
spin system to show this kind of behavior. 

We start with some simple observations. Let X be a nonexplosive conti- 
nuous-time Markov process with countable state space S (e.g., a cancellative 
spin system restricted to the space of finite states). Let X x denote the pro- 
cess X started in Xft = x and let D C S. 

Lemma 10 (Markov process leaving sets), (a) If iid xG D > s.t. X? 
D] > 0, then F[3t > s.t. X? £D] = 1 for all xeS. 

(b) // mi xeD F[3t > s.t. XI <£ D Vn > t] > 0, then F[3t > s.t. X% £ 
D Vu > t] = 1 for all x G S. 
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Proof. To prove part (b), set e := inf xeD F[3t > s.t. X% £ D Vu > t\. 
Let (J-t)t>o be the filtration generated by X x . By the Markov property and 
martingale convergence, 



(3.47) 



e!d(X s ) < F[3t > s.t. X u £D\/u>t\X s ] 
= F[3t > s.t. X u ^D\fu>t\T s \ 

► l{3t>0 s.t. X U ^D V«>t} a - s - 



This shows that lim^oo ljj(X s ) = a.s. on the complement of the event 
{3t > s.t. X u £ D \/u>t}, which implies our claim. Part (a) follows from 
part (b), applied to the process stopped at r := inf{i > : Xf £ D}. □ 

As a simple consequence of Lemma 10, we obtain the following corollary. 
Recall the definition of \VX\ in (2.10). 

Corollary 11 (Extinction versus unbounded growth). Let a < 1 and 
let X be either the neutral Neuhauser-Pacala model, affine voter model or 
rebellious voter model, and let Y be its dual parity preserving branching 
process. Then: 



(a) One has 
P[3f >0 s.t. |VXf| £{1,...,N}] = 1 

F[3t >0 s.t. \Y t y \ £{1,...,N}] = 1 

(b) If, moreover, d>2 and a > 0, then 



(3.48) 



(N > 1, \x\ < oo), 
(N > 1, \y\ < oo). 



(3.49) 



lim |VX a 



t— >oo 



oo or 3t>0 s.t. Xf 







1 



(\x\ < oo). 



Proof. We claim that, for all N > 1 and t > 0, 
(i) 



(3.50) 

and, if a > 0, 
(3.51) 



(ii) 



inf P[|VXf|^{l,...,iV}]>0, 

\Vx\<N 

inf ■ V[\Y?\£{1,...,N}]>0, 

\y\<N 



inf F\Xf = 0] > 0. 

Q<\x\<N 



Indeed, (3.50) (ii) follows from the fact that, by our assumption that a < 1, a 
particle lying sufficiently on the "outside" of y may produce N more particles 
before anything happens to the other particles of y. Formula (3.50) (i) follows 
from similar considerations. If a > 0, then by the voter model dynamics, a 
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collection of at most N ones has a uniformly positive probability to die out 
in time t, which proves (3.51). 

By Lemma 10(a), (3.50) implies (3.48). By Lemma 10(b) applied to sets 
of the form D := {x G {0, X)^ : < \x\ < N}, (3.51) implies that F[A%] = 1 
for any N > 1, \x\ < oo, where 

(3.52) A% := {3t > s.t. \X*\ > N or X? =0 Vu > t}. 
Since A X N | A x with 

(3.53) A x := ( lim \Xf\ = oo or 3t > s.t. X? = ol, 

I £ — >oo J 

we obtain Pf^P] = 1. In dimensions d > 2 this implies (3.49). □ 

In order to prove Theorem 3(a), we need to prove extinction versus un- 
bounded growth for the parity preserving branching process Y. In this case, 
Lemma 10(b) is of no use, since the analogue of (3.51) for Y does not hold 
because of parity preservation. In the following Section 3.4 we give a proof 
assuming that Y is not stable and a > 0. An alternative approach, that 
works only for very small a, but includes a = 0, is to use comparison with 
oriented percolation. For the rebellious voter model, we will use this ap- 
proach in Section 4. 

3.4. Instability. The main result of this section is the next theorem, 
which will be used in the proof of Theorem 3(a). We will also apply this 
result in [23], which is why we formulate it in some generality here. 

Theorem 12 (Extinction versus unbounded growth for parity preserving 
branching). Let Y be a spatially homogeneous parity preserving cancellative 
spin system on 7L d defined by rates a(A) satisfying (3.1). Assume that for 
L > 1, n > 0, t > 0, 



As a first step, we prove that the convergence in (3.55) holds in Cesaro mean. 

Proposition 13 (Extinction versus unbounded growth). Under the same 
assumptions as in Theorem 12, 




Assume that Y is not stable. Then 



(3.55) lim P[0 < \Y t y \ < N] = (N > 1, \y\ < oo). 



(3.56) 
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Proof. The idea of the proof is as follows. By our assumption that Y 
is not stable, one particle alone will soon produce at least three particles. In 
fact, if the process does not die out, then most of the time it will contain at 
least three particles. These particles cannot stay close together, for else they 
would annihilate each other. But single particles far from each other will 
soon each again produce at least three particles, and therefore, the number 
of particles must keep growing. 

To make this idea precise, we use induction. We write \y\ = n mod(2) to 
indicate that |y| < oo, and \y\ is even or odd depending on whether n is even 
or odd. For n > 0, consider the following statements: 



(In) 
(Hn) 

(HI n ) 



i r 1 

lim — / dtF[0 < \Y t y \ < n] = for all \y\ = n 

n ^oo T Jo 

i r T 

lim - / dt¥[3s E [0, S] s.t. < \Y t y \ < n] = 

n — »oo T Jo 

for all \y\ = n 

1 f T 

lim - / dt¥[3s £ [0, S] s.t. \Y t y +s \ =n + 2 and 

"^oo T Jo 

for some i 7^ j, \i — j\ < L] 



mod(2). 



mod(2),5>0. 



(IV, 



lim lim sup E 

S~*oo T->oo 









(0 S Jo 


T 





for all \y\ 

s 

dsl 



n 



mod(2),5>0,L> 1. 



{r n (t)<S,0<\Yy +Tn(t)+ J<n} 







mod(2), where r n (t) := inf{n >0:0<\Y t y 



for all \y\ = 

1 r T 

lim limsup — / dt- dsF[0 < \Y y \ < n] = 

5^oo T-*oo T Jo S Jo 



t+u 



<n}. 



for all \y\ 



n 



mod(2). 



We will prove that Io and Ii hold, and I n implies l n +2- Observe that if &t 
is uniformly distributed on [0,T] and independent of Y y , then 



(3.57) 



nO<\Y y \<n]=E 



dtl 



{0<\Y t y \<n} 



In the proofs below we will freely change between these and similar formulas. 

Io and Ii hold. Io is trivial. Since, by assumption, Y viewed from its lower 
left corner is not positively recurrent, the probability that Yj consists of a 
single particle tends to zero as t — > 00, which proves Ii. 

I n implies II n . This follows from the observation that 



(3.58) 



inf{P[0 < \Y t v \ < n V0 < t < 1] : < \y\ <n}=:p>0. 
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Indeed, conditional on < \ Y t y +s \ < n for some s G [0,5], with probability at 
least p, the process has between 1 and n particles during a time interval 
of length one somewhere between time t and t + 5 + 1. Therefore, if cry 
and (Ts+i are uniformly distributed on [0, T] and [0,5 + 1], respectively, 
independent of each other and of Y y , then 

limsupP[3s G [0,5] s.t. < \Y* +S \ < n] 

T->oo 

(3.59) < tttl limsupP[0 < \Y y s+1 \ < n] 

5 + 1 

■ limsupP[0 < \Yg T | < n] = 0, 



P T- 



•oo 



where in the last two steps we have used that the total variation distance 
between ¥[ax € •] and ¥[ax + crg+i G •] tends to zero as T — > oo, and our 
assumption l n , respectively. 

II n implies lll n . For n > 1, using (3.54), it is not hard to see that 

inf{P[0 < < n VI < t < 2] : \y\ = n + 2, 

(3.60) 

y(i) = 1 = y(j) for some z 7^ j, |i — j| < L} > 0. 

From this the implication follows much in the spirit of the previous impli- 
cation. This argument does not work for n = 0, so we will prove that IIIo 
holds by different means. We observe that for N = (3.54) implies that for 
L > 1 we obtain for all t > 0, 

(3.61) inf{P[F/ = 0] : y = 5i + Sj for some i / j, \i - j\ < L} > 0. 
By Lemma 10(b), this implies that, for L > 1, 



(3.62) 



lim l^y«=6i+Sj for some i^j,\i-j\<L] ~ 



1, 



which in turn implies IIIo. 

Ii, II n and lll n imply W n+ 2- Let as be uniformly distributed on [0,5], 
independent of Y y . By L and parity preservation, for each e > 0, we can 
choose So > such that, for all 5 > 5o, 

(3.63) P[|Y*|>3]>l-e. 

For each such 5, we can choose L > 1 such that 

L" 



(3.64) 



Yf°(i) = oyo<8<S,\i\ > 



> 1 - e. 



Therefore, if we start the process in a state y such that \y\ = n + 2 and 
|z — j| > L for all i 7^ j with y(z) = 1 = y(j), then with probability at least 
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(1 — 2e) n+2 , all the n + 2 particles have produced at least 3 particles at time 
05, without being influenced by each other. Thus, 



(3.65) 



F[\YV S | > 3(n + 2)] > (1 - 2e) n+2 =: 1 - e' , 



for each such y. Let a? be uniformly distributed on [0,T] and independent 
of as and Y y . If T is large, then II n and III n tell us that the probability that 



T n+2 ((TT) < S while Y y 



0"T+T n+ 2(lT T ) 



does not consist of n + 2 particles, situated 



at distance at least L from each other, is small. Therefore, by what we have 
just proved, 



(3.66) UmsupP[r n+2 ((7 r ) < S, \Y y 



T->oo 



0'T+7V 1 +2(o , t)+0'S 



< n + 2] < e'. 



Since e' can be made arbitrarily small, this proves IV n +2- 
IV n implies V n . One has 



SJo 



ds¥[0<\Y t \ s \< 



n 



E 



dsl 



Tn(t) 



{0<\Y« +3 \<n} 



(3.67) 



< E 



E 



1 rT n (t)+S 



{r n (t)<S} 



s 



dsl 



T n (t) 



{0<\Y t \J<n} 



dsl 



{r n (t)<S, 0<\Y? +Tn(t)+s \<n} 



Integrating from to T, dividing by T, and taking the limsup as T — > oo 
and then the limit S — > oo, the claim follows. 

V n implies l n . Let <tt and as be uniformly distributed on [0, T] and [0, S], 
respectively, independent of each other and of Y y . By V n , we can choose 
S(T) such that liniT^oo S(T)/T = and 



(3.68) 



IimP[0<|l^, +ff |<n]=0. 



T-^oo 



Since S(T) <C T, the total variation distance between P[ar + &s(T) € '] an d 
¥[ax G •] tends to zero as T — > oo, so I n follows. □ 

To get from Proposition 13 to Theorem 12, we need the following lemma, 
which depends on another lemma. 



Lemma 14 (Aperiodicity). Under the same assumptions as in Theo- 
rem 12, 



(3.69) 
(iV>l,5>0). 



lim sup |P[0 < \Y t y | < N] - P[0 < \Y t y +s \ < N] \ = 

t->oo ae[0,S] 
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Proof. We need to prove something like aperiodicity for an interacting 
particle system. This is not an easy problem in general. The only general 
result that we are aware of is restricted to one-dimensional systems and 
due to Mountford [20]. In our present setting, however, we can use the fact 
that each time when there are less than N particles, the next jump of our 
system happens after an exponential time with mean bounded from below, 
which causes enough uncertainty in the time variable to prove (3.69). We 
got this idea from [14], Lemma 2.4. Our way of implementing this idea is 
quite different from that reference, though. 

To make this idea rigorous, we proceed as follows. Let Z y be the embedded 
Markov chain associated with Y y , which is defined as follows. For each y £ 
{0, l} zd such that \y\ < oo, let 

(3-70) r(y) := ^ a(A)l {Ay ^ mod(2 )}. 

A 

Note that r(y) is the total rate of jumps of Y from the state y to any other 
state. It follows from (3.54) that r(y) > for all y ^ 0. Now Z y = (^) < n < noo 
is the Markov chain with state space {y 6 {0,1} Z : \y\ < oo}, started in 
Zq = y, which jumps from the state y to the state y + Ay mod(2) with 
probability a(A)/r(y). The process Z y is defined up to the random time 
rioo := inf{n > : y = 0}, which may be infinite. 

We may construct Y y from Z y as follows. Let (<r n ) n >o be i.i.d. expo- 
nentially distributed random variables with mean one, independent of Z y . 
Set 

(3.71) J7(*):=supjn>0:£a„/r(2*)<*J (i>0). 

Then rj(t) < for all t > and 

(3.72) Y t= Z n(t) (*>0)- 
Conditioning on the embedded chain Z y = (Z y )o< n < nca yields 



(3.73) 



[0 < \Y t y \ < N] = J2 P[0 < \Z y \ < N, 7](t) = n] 

n=0 

= f ¥[Z y E dz] J2 ¥ lv(t) =n\Z y = z]. 

n:0<\z n \<N 

Hence, in order to prove (3.69), it suffices to show for each z = (z n )o< n < nc 



(3.74) 



lim sup \ p [v(t) =n\Z y = z] 

~ >00 se[o,S'] n:0< ^ n | <Ar 



F[rj(t + s) = n\ Z y = z]\ = 0. 
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If < \z n \ < N for finitely many n, then 

\F[rj(t) = n\ Z y = z] - F[r)(t + s) = n \ Z y = z}\ 

n:0<\z n \<N 

(3.75) 

< p[o < z v{t) < n | z y = z] + p[o < z vit+s) <N\zy = z ], 

which tends to zero as t — > oo, uniformly for all s G [0, S]. If < \z n \ < N for 
infinitely many n, then we estimate 

\P[v(t) = n\Z y = z\- F[rj(t + s) = n \ Z y = z] | 

n: 0<\z n \<N 
oo 

(3.76) < J2 \^[v(t) =n\Z y = z]- F[rj(t + s) = n\Z y = z}\ 

n=0 

= \\F[rj(t) e-\Z y = z]- F[rj(t + s) G • | Z y = z) ||, 

where || • || denotes the total variation norm. Since < \z n \ < N for infinitely 
many n, it is easy to see that liminf n >o r(z n ) < oo. Therefore, by Lemma 15 
below, 

(3.77) lim sup \\F[r](t) G • | Z y = z] - F[r)(t + s) G • | Z y = z] \\ = 0, 
as required. □ 

Lemma 15 (Coupling of exponential variables). Let (cr n ) n >i be inde- 
pendent, exponentially distributed random variables with mean one and let 
(An)n>i be nonnegative constants. Set 

(3.78) #):=sup|n>0:f^A fcC 7 fc <t|. 
J/EfcLxA^oo, then 

(3.79) lim sup \\F[rj(t) G •] -¥[r](t + s) € -]|| =0, 

*-*°°sg[0,S] 

where || • || denotes the total variation norm. 

Proof. Our lemma is very similar to [20], Lemma 2.2, except that the 
latter uses the additional technical assumption that sup fe>1 Xk/Xk+i < oo- 
Since we do not want to assume this, our proof will be a bit different from 
the proof there. 

Let {r](t))t>o and (r]'(t))t> s be continuous-time Markov processes on N 
that jump from k — 1 to k with rate 1/A&. We start these processes at time 
and time —s in 77(0) = and rj'(-s) = 0, respectively. From time onward, 
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we let them run independently until the first time they meet, after which 
they are equal. We claim that this coupling is succesful. To see this, set 

v'(t) 

A t := ]T X k (i>0), 

k=H)(t) + l 

(3.80) r :=inf{i >0:A t = 0}, 

r R :=inf{i>0:At>ii} (R>0). 

Then (At)t>o is a nonnegative, square integrable martingale. Moreover, the 
process Mt := Af — (A)t is a martingale, where 

(3.81) (A) t := j (\( u )+i + \'(u)+i) du. 
By optional stopping, it follows that 

(3.82) E[(A) tATR ] = E[A? Ar J < E[(R V r,'(0)) 2 } (t > 0), 

so letting 1 j co, we see that (A) Tfl < oo a.s. for each R > 0. Since (A t )t>o is 
a nonnegative martingale, it has an a.s. limit as t — > oo, hence, a.s. tr = oo 
for some finite random R > 0. By (3.81), it follows that 

/"TO 

(3.83) (X v ( u )+i +\ v >(u)+i)du< oo a.s. 

Since J2kLi ^1 = °°' ^ i s n °t hard to see that X v (u)+i du = oo a.s. There- 
fore, (3.83) implies that to < oo a.s. 

This shows that our coupling is successful, hence, the total variation dis- 
tance between F[r](t) G •] and P[r]'(t) € •] tends o zero as t — > oo. It is easy to 
see that the stopping time to is stochastically increasing in s, hence, our es- 
timates are uniform for all s € [0, S], for each S > 0. Since rj(t) is distributed 
as the process in (3.78) and r}'{t) is distributed as i](t + s), our claim follows. 
□ 



Proof of Theorem 12. It suffices to show that 



(3.84) 



lim ] 

n— >oo 



[0<\Yl\<N] = 



for any sequence of times t n — > oo. By Lemma 14 and a diagonal argument, 
we can choose T n — > oo such that 

P[0 < \Y? | < N] 



(3.85) 



lim 

n— too 



1 [ n dsP[0<\Yl +s \<N] 

*-n JO 



0. 



The proof of Proposition 13 actually works more generally than for Cesaro 
times only. Let OT n and as be uniformly distributed on [0,T n ] and [0,5], 
respectively, independent of each other and of Y y . Since for each fixed S > 0, 
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the total variation distance between F[t n + ax n € •] and ¥[t n + ox n + <7s G •] 
tends to zero as n — > oo, the proof of Proposition 13 tells us that 

(3.86) n limP[0<|^ +CTr J<iV] = 0. 

Combining this with (3.85), we arrive at (3.55). □ 

3.5. Homogeneous invariant laws. In this section we prove Theorem 3, 
using Theorem 6, Corollary 11 and Theorem 12. Since Theorem 6 speaks 
about nontrivial laws, while we are interested in coexisting and nonzero laws, 
we need the following lemma. Recall the definitions of local nonsingularity 
and X-nontriviality above Theorem 6. 

Lemma 16 (Local nonsingularity). Let X be either the neutral Neuhauser- 
Pacala model, affine voter model or rebellious voter model, and let Y be its 
dual parity preserving branching process. Then, for any < a < 1: 

(a) Each homogeneous coexisting law on {0, l} zd is X -nontrivial. 

(b) Each homogeneous nonzero law on {0, 1} Z is Y -nontrivial. 

Lemma 16 can be verified by simple, but lengthy considerations, which we 
leave to the reader. To show that this point requires some care, we warn the 
reader that the claim in (a) does not hold for the disagreement voter model 
with a = 0, since for this model, the alternating configurations . . . 0101010 . . . 
are traps. 

Proof of Theorem 3. We start with part (a). Since Y survives, by 

1 /2 

Lemma 1, z/ x is not concentrated on {0,1}. Conditioning on the comple- 
ment of the set {0,1}, we obtain a homogeneous coexisting invariant law of 
X. By Lemma 16, each homogeneous coexisting invariant law of X is X- 
nontrivial in the sense of Theorem 6. Therefore, since a < 1, it follows from 

1/2 

Theorem 6 and Corollary 11(a) that ^ x is the only homogeneous coexisting 
invariant law. If, moreover a > and Y is not stable, then by Theorem 6 
and Theorem 12, the convergence in (1.13) holds. Here, condition (3.54) of 
Theorem 12 is easily seen to follow from our assumption that a > 0. 

To prove also part (b), we observe that since X survives, by Lemma 1, 
Vy 2 is not concentrated on 0, hence, Vy 2 conditioned on being nonzero is a 
homogeneous nonzero invariant law of Y. By Lemma 16, each homogeneous 
nonzero invariant law of Y is Y-nontrivial. Therefore, since a < 1, it follows 

1/2 

from Theorem 6 and Corollary 11(a) that Vy is the only homogeneous 
nonzero invariant law. If, moreover, a > and d>2, then by Theorem 6 
and Corollary 11(b), the convergence in (1.14) holds. □ 
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3.6. A counterexample to a result by Simonelli. Our Theorem 6 is similar 
to Theorem 1 in [22]. The proof is also similar, with his Theorem 2 playing 
the same role as our Corollary 9. An important difference is that while we 
use the "norm" || • ||g defined in (3.11), Simonelli works with the usual l 1 - 
norm | • |. As we will see in a moment, the result of this is that his Theorem 2 
is false. We do not know if his main result Theorem 1 is correct or not; as 
it is, his proof depends on his Theorem 2, and is therefore not correct. 

More precisely, in [22] it is assumed that X and Y are interacting particle 
systems on Z rf satisfying a duality relation of the form (1.8). It is assumed 
that 



where Y y denotes the process y started in Yq = y. Then, translated into our 
terminology, Simonelli ([22], Theorem 2) states that: 

If P[Yo £ •] is homogeneous and nonzero, and x n £ {0, 1} Z satisfy linin^,^ |a;„| = 
oo, then lim„_ K x>P[|i / ta;n| is odd] = ^ for all t > 0. 

This claim is false, as is shown by the following: 

Counterexample. We pick some a € (0, 1) and we take for Y the 
ADBARW, which is the dual of the rebellious voter model. It is easy to see 
that Y satisfies (3.87). Set P[io = i] = 1> which is homogeneous and nonzero, 
and x n (i) := 1 for i = 1, . . . , n and x n (i) := otherwise. Using duality, we see 
that P[|YfX n | is odd] = P[|l rj Af n | is odd], where X Xn denotes the rebellious 
voter model started in Xq" = x n . The dynamics of the rebellious voter model 
are such that at t = 0, the only places where sites can flip are at the endpoints 
of the interval {1, ...,n}. Therefore, it is easy to see that ¥[X^ n = x n ] > 
1 — e _4i , and as a result, 



where the convergence is uniform in n. It follows that, for t sufficiently small, 
the limit linv^oo P[|Y<x n | is odd] does not exist. 

There seems to be no easy way to repair Simonelli's Theorem 2. The 
essential observation behind our Lemma 7 is that, instead of x n being large, 
one needs that x n contains many places where parity can change. For the 
rebellious voter model, this means that x n must contain many places where 
the two types meet. Our proof of Lemma 7 differs substantially from the 
methods used in [22]. 

The place where Simonelli's proof goes wrong is his formula (5), where it 
is claimed that if P[io G "] is homogeneous and nonzero, then for each e > 
and t > there exists a 5 > such that, for all s € [0,e], 



(3.87) 



0<P[y^(i) = l]<l (t>0,i£Z d ,y^0) 



(3.88) 




if n is odd, 
if n is even 



F[Y t €{y:5< F[Y s y {0) = 1] < 1 - 5}] > 1 - e 



([22], (5)). 
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While this inequality is true for any s > fixed with 5 = 8(t, s, e) depending 
on s,t, and e due to his assumption (3.87) and continuity of probability 
measures, 5 cannot be chosen uniformly in s£ [0,e]. In fact, for s = 0, the 
set {y : 5 < P[17(0) = 1] < 1 - 6} is empty for any 5 > 0. 

4. Complete convergence of the rebellious voter model. In this sec- 
tion we prove coexistence and complete convergence for the rebellious voter 
model, as stated in Theorem 4. The main tool for this will be Theorem 5, 
which states that the dual ADBARW dominates oriented percolation. This 
theorem is formulated and proven in the following Section 4.1. 

4.1. Comparison with oriented percolation. We let Y be an ADBARW 
started in an arbitrary deterministic initial state Yq = y € {0, 1} Z , and we 
define sets of "good" points (xn)n>o as in (2.13). We start by considering 
a single time step in the case that a = 0. Our first result says that "good" 
events have a large probability. 

Proposition 17 (Good events are probable). Assume that a = 0. Then, 
for each p < 1, there exist L>1 and T > 0, not depending on the initial state 
y, such that F[x € xi] — P for all x € Z Q dd such that xo H — 1, x + 1} ^ 0. 

Proof. Our basic observation is that, in case a = 0, the right-most par- 
ticle of an ADBARW started in a finite initial state has a drift to the right. 
Indeed, if rt '■= max{i € Z: Yj(i) = 1}, then depending on the configuration 
near the right-most particle, the changes in rt due to the various possible 
jumps and the resulting drift are as follows: 

configuration change in rt 
...00100... +2 

...01100... 
...10100... 

...11100... 

This shows that in each configuration the drift is at least one. We note, 
however, that it is not possible to stochastically bound rt from the left by 
a random walk with positive drift (independent of anything else). This will 
cause a slight complication in what follows; in fact, we will use two random 
walks that become active when Y(rt — 1) = or 1, respectively [see formula 
(4.11) below]. 

We need to prove that if xo H {x — 1, x + 1} ^ for some x € Z D dd, then 
¥[x £ xi] ^ P- By symmetry, we may without loss of generality assume that 



rate drift 


1 


+2 


1] 










1] 










!! 
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x = 1 and S xo- To simplify notation, let us identify subsets of Z with 
their indicator functions. Then, assuming that y n {— L, . . . , L} ^ 0, which 
is equivalent to G Xo> we need to show that the probability 

(4.2) F[Y T n {L, . . . , 3L} + and Y t n {-2L, . . . , 6L} / VO < t < T] 

can be made arbitrarily large by choosing L and T appropriately. In view 
of this, we are actually not interested in the right-most particle of our AD- 
BARW, but in the particle that is closest to our target 2L. Thus, we put 

(4.3) s t := inf{i > : Y t (2L - i) V Y t (2L + i) = 1}. 

Assuming that so < 3L which follows from y n {— L, ... ,L} ^ 0, we need to 
show that the probability 



(4.4) 



st < L, sup st < 4L 

0<t<T 



can be made arbitrarily large. For any n > 0, we set 

r<„ :=inf{t >0:s t < n}, 

(4.5) 

r>„ := inf{t > : S( > n}. 

We observe from (4.1) that whenever an ADBARW borders at least two 
empty sites, it tends to invade these with a drift of at least one. In view of 
this, we choose T = 2L. By Lemmas 18 and 19 below, there exist constants 
C, A > such that 

(4.6) P[r<2 < T and t> 4 l > T] > 1 - Ce~ Ai . 

Using Lemma 19 once more, we see that, moreover, for some C", A' > 0, 

(4.7) F[s t < L Vr< 2 < t < T] > 1 - C'e~* L . 

Combining these two estimates, we see that the probability in (4.4) can be 
made as close to one as one wishes by choosing L large enough. □ 

We still need to prove two lemmas. 

Lemma 18 (Hitting the target). For each 5 > 0, there exist constants 
C, A > such that if so < K , then 

(4.8) P[r< 2 > (1 + 5)K] < Ce~ XK (K > 1) . 

Proof. Let 

(4.9) :=inf{i>0:Y t (2L-i) = l}. 
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Since the evolution rules of the process are symmetric, we may assume with- 
out loss of generality that Sq < K. We set 

(4.10) <j> q {t):= f V(.-_i)= g} du (t>0,g = 0,l). 

Let RP^R 1 be continuous-time random walks on Z, starting in zero, with 
the following jump rates: 

random walk jump size rate 

(4.11) 

Rl 

In view of (4.1), we can couple sj to R?,R} in such a way that 

(4.12) st < so - R° Mt) - R\ l{t) (0 < t < r< 2 ), 

where R° and R 1 are independent of each other and of <fio, (f>\. It follows from 
large deviation theory, more precisely, from Cramer's theorem (Theorem 27.3 
in [16]) and a little calculation, that for each e > there exist constants C £ 
and X £ > such that 

(4.13) F[\R?-t\ >et] <C £ e- Xst (t>0,g = 0,l). 

I Ag su 

(4.14) ¥[\R q s - s\>st]<C' e e- K - t (0 < s < t, q = 0, 1). 



+2 


1 


-1 


1 


+2 


1 


+1 


1 


-2 


1 



We claim that there exist C' e and X' £ such that 



?1 



Since it is not obvious that P[|i2f — s\ >et] < P[|i?j — 1\ > et], it is not entirely 
trivial to get from (4.13) to (4.14). Here is a clumsy argument: If < s < t, 
then 

(4.15) F[\R q s -s\> et] < F[\Rj -s\> es] < C £ e-^ s < C £ e- {1/2)x ^. 

If < s < ^t, then, by the independence of random walk increments and 
what we have just proved, 

F[\R q t -t\> \et] > F[\Rj -s\> et]P[\R?_ s -(t-s)\< \et] 

(4.16) 

>P[|^- S |>rf](l-C £/2 e-( 1 /2)A s/2 *) ) 

hence, by (4.13), 

(4.17) F[\Rj -s\> et] < (1 - C^e-^ 2 ^)- 1 ^-^*. 
Combining (4.15) and (4.17), we obtain (4.14). 
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To prove (4.8), we set M := (1 + 6)K, we choose e such that (l-e)M = K, 
and observe that by (4.12), the fact that </> (M) + 0i(M) =M, and (4.14), 

P[r< 2 > M] < F[R° M m) + R MM) < (1 " e)M] 

(4-18) < n\R° M M) + 4i(M) " MM) - MM)\ > eM] 

<J2n\R q MM ) - MM)\ > \sM) < 2C' £/2 e- X '^ M . 
Setting C := 2C' e/2 and A := X' e/2 (1 + 6), we arrive at (4.8). □ 

Lemma 19 (Escaping the target). For each 5> 0, there exist constants 
C, A > such that if sq < K, then 

(4.19) nr>{i+5)K<T\<C{T+l)e~ XK (T>0,K>1). 

Proof. We start by showing that we can choose A > such that if so < 
L' <L, then 

(4.20) P[r< 2 >r> L ]<e- A(L - i ' ) . 

Indeed, it is not hard to see from (4.1) that, for A sufficiently small, the 
process 

(4.21) M t := e A ' inf{ *' T <2 } (t > 0) 

is a super martingale. Setting r := inf{r<2, t>^}, by optional stopping, it 
follows that 

(4.22) e AL P[s T > L] < E[M T ] < e XL ' , 

which proves (4.20). To prove (4.19), we note that each time the process 
St enters {0,1,2}, it must stay there at least an exponentially distributed 
time with mean i. Therefore, the number of excursions from {0, 1, 2} during 
a time interval of length T is bounded by a Poisson random variable with 
mean 2T, and by (4.20), the number of excursions from {0,1,2} that go 
beyond (1 + 5)K is bounded by a Poisson random variable W with 

(4.23) P[W > 0] < E[W] = 2Te~ m+5)K - 2) . 

Using (4.19) once more, we conclude that the probability in (4.19) is bounded 

b y e -X5K +2Te -X((l+5)K-2)_ n 

We now turn to the proof of Theorem 5. It suffices to prove the statement 
for a single time step; the general statement then follows by induction. Thus, 
we need to show that for each p < 1 we can find L,T so that we can de- 
fine i.i.d. Bernoulli random variables {^(^,1) : % G ^odd} with Pfav^i) = l]=p 
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such that x G xi whenever u x i = 1 and xo H {x — l,x + 1} ^ 0. Since we 
can thin the set of points where U( x ,i) = 1 if necessary, it suffices if the 
{u( x ,i) -x G Z dd} are independent and F[uj( x ^ = l]>p for each x. It actu- 
ally suffices if the '■ x € %odd} are m-dependent for some fixed m > 1, 
since a well-known result (see [19], Theorem B26) tells us that m-dependent 
random variables with intensity p can be estimated from below by indepen- 
dent random variables with intensity p' =p'(p,m) depending on p and m 
in such a way that lim p ^i p'(p,m) = 1. See also [12] who considers a more 
general form of m-dependence for oriented percolation and the comparison 
argument. 

Proof of Theorem 5. As mentioned before, it suffices to prove the 
statement for a single time step, and for the process started in any deter- 
ministic initial state Yq = y. A naive approach is to put W(z,i) '■= ^{xexi} ^ or 
x G J := {x G Z: xo H {x — l,x + 1} ^ 0} and U( x ,i) := 1 otherwise. However, 
since these events are not m-dependent for any fixed m, we are going to 
extend our definition of a "good" event, in such a way that the new events 
still have a high probability, and are m-dependent. To that aim, we put 

Xi '■= {x G xi : there is no path in the graphical representation ending 

(4.24) at time T in {2Lx - 4L, . . . , 2Lx + 4L} and starting 

at time outside {2Lx — 11L, . . . , 2Lx + 11L}}, 

where, as before, we choose T = 2L. The motivation for this is as follows. 
Let 

r t := sup{i S Z : there is a path starting at (i, T — t) 

(4.25) 

and ending in {T} x {2Lx - 4L, . . . , 2Lx + 4L}}. 

It is not hard to see that rt can be bounded from above by a random 
walk that makes jumps of size +1 and +2, both with rate 1. Therefore, 
the expected distance covered by such a random walk is 3T = 6L. A large 
deviation estimate of the same sort as used in the proof of Lemma 18 now 
tells us that for each x € Z oc jd , the probability of there being a path ending 
at time T in {2Lx — 4L, . . . , 2Lx + 4L} and starting at time outside {2Lx — 
11L, . . . , 2Lx + 11L} tends to zero exponentially fast as L — ► 00. Combining 
this with Proposition 17, we see that for each p we can choose L,T in such 
a way that F[x € x'i] — P f° r all x £ J. Choosing a' close enough to zero so 
that the probability of any event with rate a' happening in the graphical 
representation in the block [0, T] x {2Lx — 11L, . . . , 2Lx + 11L} is small, we 
conclude that 

For each p < 1 there exist an a' > such that for all a G 
(4.26) [0, a') there exist L > 1 and T > such that F[x G x'i] > P 
for all x G J. 
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We now put ui( x ,i) '■= ^{xex[} ^ or x ^ J and <^( x ,i) := 1 otherwise. The {^(a;,i) : x £ 
^odd} constructed in this way are m-dependent for a suitable m [in fact, 
^( x ,i) an d are independent if \x — x'\ > 12], so by the arguments pre- 

ceding this proof, they can be bounded from below by i.i.d. Bernoulli random 
variables with an intensity that can be made arbitrarily high. □ 



4.2. Complete convergence. In this section we prove Theorem 4 about 
coexistence and complete convergence for the rebellious voter model X for 
small enough a. Throughout this section we choose some pG (p c , 1), where 
p c is the critical value for survival of oriented percolation, and we fix a' > 0, 
L > 1, and T > such that the process (xn)n>o defined in (2.13) can be 
coupled to an oriented percolation process (W n ) n >o with parameter p as in 
Theorem 5. The proof of Theorem 4 is based on the following two lemmas. 
Recall from Section 2.1 that the ADBARW is both the dual and the interface 
model of the rebellious voter model. 



Lemma 20 (Almost sure extinction versus unbounded growth). For all 
a £ [0, a'), the ADBARW Y y started in any finite initial state Yq = y satis- 
fies 



(4.27) 



lim |y/| 

t— »oo 



oo or 3t > s.t Y t v = 



1 



(\y\ <oo). 



Moreover, Pflim^oo \Y^\ = oo] > for all y ^ 0. 



Lemma 21 (Intersection of independent processes). For all a € [0, a'), 
two independent ADBARW s Y and Y' satisfy for t > 0, N > 1 

(4.28) lim F[\Y t Y/\ >N}= F[Y t + Vi > 0] • P[Y/ ^ Vt > 0]. 

t^oo 

Proof of Theorem 4. If a is as in Lemma 20, then the ADBARW is 
unstable and survives. Therefore, by Lemmas 1 and 2, the rebellious voter 
model exhibits coexistence, survival, and its dual ADBARW is not stable. 
This proves part (a). 

To prove part (b), it suffices to consider deterministic initial states. Let 
X x denote the rebellious voter model started in Xq = x and set p q (x) := 
¥[Xf = q for some t > 0] (q = 0, 1). By (2.4), it suffices to show that, for 
each \y\ < oo, 

lim F[\Xfy\ is odd ] 

(4.29) 

= (1 - p (x) - Pl {x))\T[Y^ V< > 0] + P i(x)l M is odd} . 
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Since the dynamics of X are symmetric for the O's and l's, we may without 
loss of generality assume that the starting configuration x has infinitely 
many O's. In this case p\{x) =0 and (4.29) reduces to 

(4.30) lim F[\Xfy\ is odd ] = i(l - p {x))¥[Y t y ^ Vt > 0]. 

t— >oo 

Since, by duality, 

F[\Xfy\ is odd ] 

(4.31) 

' P[*ft-i)/ 2 G _/ P[^-i)/ 2 e dy]P[\XfYy\ is odd], 
by Lemma 7, in order to prove (4.30), it suffices to show that 
nmP[|{ieZ:rJ_ 1)/2 (i)=l, 

(4.32) *Tt-i)/ 2 (* + 1) + *(t-i)/ 3 (* + 2)}| > AT] 

= p[xf ^ o vt > o] • p[y/ / o vt > o] (iv > i). 

Let y/(i) := 1{ x x (i+i)^x x (i+2)}) that is, F' is the interface model of X x , 
translated over a distance one. Then (4.32) simplifies for N > 1 to 

(4.33) Hm P[|^_ 1)/2 yJ_ 1)/2 1 > N] = P[F/ ^ Vt > 0] • F[Y t y ± V* > 0], 
which is true by Lemma 21. □ 

The rest of this section is occupied by the proofs of Lemmas 20 and 21. 

Proof of Lemma 20. By Lemma 10(b), the claims will follow provided 
that we show that 



(4.34) inf P 

|»l>o 



lim \Yt\ = 00 

t— >oo 



>0. 



By translation invariance, it suffices to consider the infimum over all y such 
that y(0) = 1. Choose p> p c , the critical value for survival of oriented perco- 
lation. By Theorem 5, the process (xn)n>o defined in (2.13) can be coupled to 
an oriented percolation process (W n ) n >0 such that Wo = xo and W n C Xn for 
all n > 1. Since G xo due to the fact that y(Q) = 1, and since p> p c , there is 
a positive probability that lim ra _» 00 \W n \ = oo and, hence, limt^oo \Y y \ = oo. 
□ 



The proof of Lemma 21 is somewhat more involved. We start with some 
preparatory lemmas. Our first lemma says that if y and y' are close in many 
places, then |1^V| is large with probability close to one. 
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Lemma 22 (Charging target sets). Assume that a < 1, let Y y be an AD- 
BARW started in y, let y' G {0, 1} Z , and let K > 1, t > 0. Set D K (y,y') ■= 
G 1? :y(i) = 1 = y'(j), \i-j\<K}. Then 

(4.35) lim inf F[\Y t v y' \ > N] = 1 (N>1). 

M->oo\D K (y,y')\>M 

Proof. Set C := {j:3i s.t. G D K (y,y')}. For each j G C, choose 
in some unique way a site i with y(i) = 1 for which \i — j\ is minimal. Let 
I := {i, . . . ,j} if i < j and / := {j, . . . , i} if j < i. Let Gj denote the event 
that in the graphical representation for Y y (see Section 3.1), there is an odd 
number of paths from (i,0) to (j,t), while during the time interval [0,t], 
there are no arrows starting outside / and ending in /. Then, for given K 
and t, the probability of Gj is uniformly bounded from below. To see this, 
by symmetry, we may assume i < j. Then the particle at i may branch to 
the right, producing two particles at i + 1 and i + 2, which can again branch 
to the right, creating, in a finite number of steps, a particle at j, while with 
positive probability, nothing else happens in /. Now, if \DK(y,y')\ > M, 
then we can select C'cC such that \j — j'\ > 2K + 1 for each G C 
with j / f, and \C'\ > M/(2K + l) 2 . Then the events Gj with j G C are 
independent with a probability that is uniformly bounded from below, hence, 
if M is sufficiently large, then with large probability many of these events 
will occur. This proves (4.35). □ 

In what follows, for any x G Z, we define I x as in (2.12), and for y G {0, 1} Z , 
we define 

(4.36) r){y) := {x G Z cvcn : 3i G 4 s.t. y(t) = 1}. 



Lemma 23 (Charging target intervals). If Y is an ADBARW with pa- 
rameter a G [0,a'), started in Yq = y, then 

(4.37) lim F[Y t7 t and \ V (Y t )\ < N] = (N>1). 

t— >oo 

Proof. By Lemma 22 applied with K = L and y' = J2x<=z $2Lx, for each 
t > 0, 

(4.38) lim inf P[|7?(y/)| > iV] = 1 (N > 1). 

M— >oo |j/|>A/ 

Now write 

F[Y t ^0 and |»/(y t )| < TV] 

(4.39) =P[Y t t^O and |r?(y t )| < TV | < |Y t _i| < M]P[0 < |Y_i| < M] 

+ P[Y^0 and |T/(r t )|<iV| |l^_i|>M]P[|n_i|>M]. 
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By Lemma 20, the first term on the right-hand side of (4.39) tends to zero 
as t — > oo, while by (4.38), the limsup as t — > oo of the second term can be 
made arbitrarily small by choosing M large enough. □ 

Fix Bernoulli random variables {u> z : z € 2^ vcn } with intensity p as in Sec- 
tion 4.1, and for each A C Z even , let W A = (W£) n >o denote the oriented 
percolation process started in A defined in (2.11). Using the same Bernoulli 
random variables u z , we can define a process W = (W n ) n& z by 

(4.40) W n :={xeZ:(x,n)eZ 2 cvcn ,-oo^(x,n)} (n€Z), 

where — oo — » (x,n) means that there exists an infinite open path with re- 
spect to the co z , starting at time — oo and ending at (x,n). Then W is a 
stationary (with respect to shifts on Zg ven ) oriented percolation process. We 
call 

(4.41) %:=P[F 2 n£-] (n€Z) 

the upper vfivQ/ricLTtt law of W '. It is known tha/t, for each 7i ^ 1 and. x £ ^ even , 

(4.42) ItoP^ + and W n n [-K, K] £ W^ x} ] = 0. 

This follows, for example, from [17], Theorem 2.27 (see also Theorem 2.28). 
The statements there are for the one-dimensional nearest neighbor voter 
model, but the proofs apply to our setting as well. Alternatively, one may 
consult [10] and [2], Theorem 5, where this (as well as the much more power- 
ful shape theorem) is proved in a multidimensional setting, for p sufficiently 
close to one and for p> p c , respectively. 

Lemma 24 (Local comparison with upper invariant law). Let Y be an 
ADBARW with parameter a G [0,c/), started in Yq = y, and let t^ — > oo. 
Then, for each k, we can couple Y tk to a random variable with law z7 W; 
in such a way that 

(4.43) limmfF[V k n[-K,K]Cri(Y tk )]>F[Y t ^0 Vt > 0] (K>1). 

k— »oo 

Proof. By Theorem 5, for each s > 0, we can couple Y to an ori- 
ented percolation process W s = (W^) n >o started in Wq = n(Y s ), defined 
by Bernoulli random variables that are independent of (Y u ) ue [ 0>s i, in such a 
way that 

(4.44) W°C?i(Y s+nT ) (n>0). 

We couple W s to an "upper invariant" oriented percolation process W s as 
in (4.40), and put V SyTl := W\ n (n > 0). Fix K > 1 and e > 0. It is not hard to 
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see that lrni/v-»oo hrflA^ArPp^f ^ Vn > 0] = 1. Therefore, by Lemma 23, 
we can choose so > such that, for all s > sq, 

(4.45) P[W* + Vn > 0] > F[Y S ^ 0] - \e. 

Since the process Y cannot move infinitely far in finite time, it is not hard 
to see from (4.42) that for each s£ [sq, sq + 2nT), we can choose no > such 
that, for all n > no, 

(4.46) F[V s , n n [-K, K] C n(y s+2nT )] > P[I^ n s + Vn > 0] - \e. 

It follows that for all t > sq + 2noT we can choose s S [so> s o + 2nT) and 
n > no such that s + 2nT = i and 

(4.47) F[V s>n n [-if, if] c r?(y t )] > PpS ^ Vt > 0] - e. 

Since X and e are arbitrary, for each — > oo, we can find — ► oo, — > 0, 
and Sk,rik such that 

(4.48) p[v Sfc>nfe n [-if^ift] c r?(y tfc )] > prr t ^ o vt > o] - £fc , 

which proves our claim. □ 

Proof of Lemma 21. Set p := F[Y t / \ft > 0] and pt := F[Y{ f Vt > 
0]. It suffices to prove the convergence in (4.28) along an arbitrary sequence 
of times — > oo. It is clear that the limsup is bounded from above by pp' . 
To bound the liminf from below, by Lemma 22, it suffices to prove that 

(4.49) limMF[D 2L (Y tk ^,Y')>N]>p P ' (N>1). 

k— >oo K 

By Lemma 24, we can couple (Y tk -\, YJ) to random variables Vk,V£ each 
having law I7 W , in such a way that (Y tk -i,Vk) is independent of (Y t ' ,V£), 
and for if > 1, 

(4.50) lhninf F[V k n [-if, if] C r^-i),^ D [-if, if] C r?(Y t J] > pp'. 

It is not hard to see that if V, V is a pair of independent random variables 
each having law i/\y, then \V PI V'| = oo a.s. Therefore, for each N, e > 0, we 
can choose if large enough such that 

(4.51) F[\VnV n[-K,K]\<N]<e, 
and hence, by (4.50), 

(4.52) \MF[D 2L {Y tk ^X k )>N}>pp' -e. 
Since e > is arbitrary, (4.49) follows. □ 
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